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Abstract 

In this paper the percolation of monomers on a square lattice is studied as the 
particles interact with either repulsive or attractive energies. By means of a finite- 
size scaling analysis, the critical exponents and the scaling collapsing of the fraction 
of percolating lattice are found. A phase diagram separating a percolating from a 
non-percolating region is determined. The main features of the phase diagram are 
discussed in terms of simple considerations related to the interactions present in the 
problem. The influence of the phase transitions occurring in the system is reflected 
by the phase diagram. In addition, a scaling treatment maintaining constant the 
surface coverage and varying the temperature of the system is performed. In all the 
considered cases, the universality class of the model is found to be the same as for 
the random percolation model. 
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1 Introduction 



Percolation theory has attracted a great deal of interest in the last few decades and the 



activity in the field is still growing PElEllllElEllTllHlIHlIIOllIIlliaiTS]. This is mainly 



because some aspects of the percolation process such as the geometrical phase transitions 
occurring in the system have gained a particular impetus due to the introduction of 
techniques such as Monte Carlo (MC) simulations and series expansions [HKIlj. However, 
the problem is far from being exhausted. 

In fact, most of the studies are devoted to the percolation of molecules that are ir- 
reversibly deposited. In part this is due to the fact that the deposition (or irreversible 
adsorption) of particles on solid surfaces is a subject of considerable practical importance. 
In many experiments on adhesion of colloidal particles and proteins on solid substrates, the 
relaxation time scales are much longer than the times of the formation of the deposit. In 
such processes, the temperature of the system does not play any relevant role and it is not 
considered. However, in numerous systems of both theoretical and practical importance, 
where the adsorbed particles are in thermodynamic equilibrium, the spatial distribution 
of the adsorbate might be characterized by using the percolation model fibi [16] . In these 
cases, the temperature governs the phase in the system and can be an important con- 
trolling factor in the percolation process. In the simplest case, where repulsive monomers 
are adsorbed on a square lattice, the system exhibits a continuous phase transition from 
a disordered state at high T to a doubly-degenerated c(2 x 2) ordered state at low T, 
with a critical temperature which satisfies Kc = 1.7637 [17], with K = w/kBT, ks the 
Boltzmann constant and w the interaction energy between nearest-neighbour adatoms. 
On the other hand, the same system with attractively interacting particles goes through 
a first-order phase transition. It is clear that a percolation study on the spatial configura- 
tion of the adlayer should recognize peculiarities related to the phase transitions occurring 
in the adsorbate. Such a study implies the determination of the critical parameters as 
a function of the concentration (surface coverage) and temperature. To the best of our 
knowledge this study has not been done and it is the main purpose of this paper. 

The paper is organized as follows. In section [2] the model of adsorption of interact- 
ing monomers on a two-dimensional square lattice is presented. The analysis of results 
obtained by using finite-size scaling theory is given in section [31 In section [H, the phase 
diagram is discussed along with the basis of a thermal finite-size scaling study. Finally, 
conclusions are drawn in section [5] 

2 Basic definitions 

Let us consider that the substrate is represented by a two-dimensional square lattice of 
M = L X L equivalent adsorption sites, with periodic boundary conditions. In order to 
describe the system of N monomers adsorbed on M sites at a given temperature T, let 
us introduce the occupation variable q which can take the following values: 



* 1 0, if site i is vacant. 

Particles can be adsorbed on the substrate with the restriction of at most one adsorbed 
particle per site and we consider a nearest-neighbour (NN) interaction energy w between 
them. Under these considerations, the Hamiltonian of the system is given by 
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M 



H = w 



(2) 



where represents pairs of A^A^ sites and Si is the adsorption energy of the sites on 

the surface. In addition, we have taken = e = without loss of generahty. 

For fixed values of surface coverage, 9 = N/M, and temperature T, the thermody- 
namic equilibrium is reached in the canonical ensemble by using a standard Kawasaki 
algorithm [18]. The procedure is as follows. An initial arbitrary configuration of A^ ad- 
sorbed monomers with the desired surface coverage is generated. Two sites are randomly 
selected, and their positions are established. Then, an attempt is made to interchange 
their occupancy state with probability given by the Metropolis rule [I9j : 



where AH = Hf—Hi is the difference between the Hamiltonian evaluated at the final state 
and the one computed at the initial state. A Monte Carlo step (MCS) is achieved when M 
pair of sites have been tested to change its occupancy state. The time (expressed in units 
of MCS) required for equilibration depends on lattice size, temperature and coverage. 
Typically, m = 10^ MCSs suffices for smaller lattices containing up to 32 x 32 sites in the 
whole range of both temperature and coverage. Then, a set of m = 2 x 10'^ samples in 
thermal equilibrium is generated by taking configurations separated from each other by r 
MCSs in order to avoid possible correlations between the states. In the smallest lattices 
considered in this paper, r = 10^ prevents such undesired correlations. The accuracy of 
the procedure and the correctness of the algorithm were tested by obtaining the behaviour 
of different quantities (for example, the adsorption isotherms (surface coverage, 6, versus 
the normalized chemical potential n/ksT), the specific heat, the order parameter, etc) 
and comparing with the corresponding ones derived from the real space renormalization 
group approach (RSRG) [20]. 

The central idea of the percolation theory is based on finding the minimum concen- 
tration 9 for which at least a cluster (a group of occupied sites in such a way that each 
site has at least one occupied nearest-neighbour site) extends from one side to the oppo- 
site one of the system. This particular value of the concentration rate is named critical 
concentration or percolation threshold and determines a phase transition in the system. 
In this paper, the percolation process will be studied under two different perspectives. 
Namely, samples will be prepared for fixed temperature (coverage) and variable coverage 
(temperature). We call this feature coverage percolation {thermal percolation) . 

In the random percolation model, a single site is occupied with probability p. The 
samples are then generated through irreversible adsorption. In our problem, the occu- 
pancy state of every site strongly depends on K. In both cases, for a precise value of 
concentration, the percolation threshold of sites, at least one spanning cluster connects 
the borders of the system (indeed, there exists a finite probability of finding n {> 1) 
spanning clusters [2T1 [221 E21 [21] ) • Then, a second-order phase transition appears at such 
coverage which is characterized by well-defined critical exponents. It should be empha- 
sized that for K = the system becomes uncorrelated and formally we have the random 
percolation model. 

A study of the finite-size effects allows us to make a reliable extrapolation to the 
thermodynamic limit (L oo). Details of this study will be given below. 



P = min {1, exp {—AH/kBT)} 



(3) 
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3 Finite-size scaling 



It is well known that it is a quite difficult matter to analytically determine the value of the 
percolation threshold for a given lattice [21 |5l El El E] • For some special types of lattices, 
geometrical considerations enable us to derive their percolation thresholds exactly. Thus, 
exact thresholds for the random percolation problem are known for (a) square, triangular 
and honeycomb lattices and (b) triangular and Kagome lattice concerning the bond and 
site problem, respectively. In both cases, analytical results are obtained when a monomeric 
species is considered. For different conditions, i.e. for systems which do not present such a 
topological advantage, percolation thresholds have to be estimated numerically by means 
of computer simulations. 

As the scaling theory predicts [nj, the larger the system size to study, the more accu- 
rate the values of the threshold obtained therefrom. Thus, the finite-size scaling theory 
gives us the basis to achieve the percolation threshold and the critical exponents of a sys- 
tem with a reasonable accuracy. For this purpose, the probability R = Ri{0) that a lattice 
composed of L x L elements (sites or bonds) percolates at concentration 9 can be defined 
[2]. Here, as in [251 [26], the following definitions can be given according to the meaning of 
X: (a) R^^^\6) is the probability of finding a rightward (downward) percolating cluster; 
(b) Ri{9) is the probability that we find a cluster which percolates both in a rightward 
and in a downward direction; (c) R^{0) is the probability of finding either a rightward or a 
downward percolating cluster and (d) Rt{e) = \ (^) + i?f (^)] = \ [i?i(^) + i?^(^)] . 

The ffist step for determining the percolation threshold consists in evaluating the 
effective threshold QJ^V) (the concentration at which the slope of Rl{0) is the largest) for 
a lattice of finite-size L. In the MC simulations, Rf^{0) is determined for each discrete 
value of 6 according to the considered finite regular lattice [13]. Once the temperature 
is fixed, the next procedure is followed: (a) the construction of m samples for a given 
coverage (according to the scheme presented in section [2]) and (b) the cluster analysis by 
using the Hoshen and Kopelman algorithm [27]. In the last step, the number of clusters 
for each sample, n^, of size s (a cluster of size s is composed by s connected elements) is 
determined in order to verify whether a percolating island exists. This spanning cluster 
could be determined by using the criteria R, D, I ot U. m runs of two such steps are 
carried out for obtaining the number of them for which a percolating cluster of the 
desired criterion X is found. Then, Rl{0) = /n is defined and the procedure is 
repeated for different values both of 9 and lattice sizes, L {L = 32, 48, 64, 80, 96 and 128). 

In figure 1, the probabilities R^iO) (squares), R^{0) (circles) and R^iO) (triangles) are 
presented. Three different values of K are shown. From a first inspection of the figure 
(and from data not shown here for the sake of clarity) it is observed that (a) curves cross 
each other in a unique universal point, R^ , which depends on the criterion X used; (b) 
those points do not modify their height for the different Ks. This finding indicates, as 
is expected, that the universality class of the phase transition involved in the problem is 
conserved no matter what the value of K is; (c) those points are located at very well- 
defined values in the 6'-axes, determining the critical percolation threshold for each K\ 
(d) the standard percolation problem is recovered for = giving a critical coverage 
Oc = 0.5927 and (e) 9c increases (decreases) for positive (negative) value of K. A detailed 
explanation of this point will be given below. 

The second step is the extrapolation of 0^{L) towards the limit L ^ oo by using the 
scaling hypothesis. Thus, the correlation length, ^ , can be expressed as 
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(4) 



where the critical exponent v is analytically shown to be equal to z/ = 4/3 in the case of 
random percolation [21 IH El E]- As = 0^{L) the correlation length reaches the linear 
dimension L of the lattice. Thus, we have 



where A is a non- universal constant. Figure 2 shows the extrapolation towards the 
thermodynamic limit of 9^{L) according to equation ([5]) for different values of K as 
indicated. This figure lends support to the assertion given by equation ([5]): (a) all the 
curves (different criteria) are well correlated by a linear function, (b) they have a quite 
similar value for the ordinate in the L — > oo and (c) the fitting determines a different value 
of the constant A depending on the type of criterion used. It is also important to note that 
0^{L) gives an almost perfect horizontal line which is a great advantage of the method 
because it does not require precise values of critical exponents in the process of estimating 
percolation thresholds. The maximum of the differences between \9^^{oo) — 6*^(00)! and 
I^^(cxd) — 9^{oo)\ give the error bar for each determination of 9^. 

The scaling law hypothesis also predicts the collapsing of the curves R^{9) when they 
are plotted as a function of a reduced variable u = {9 — 9c) L^^'^: 



R-^{u) being with the scaling function. Thus, is a universal function with respect to 
the variable u. In figure 3(a), as an illustration, we plot R^ as a function of u for K = 0. 
This gives an additional indication for the numerical value of the critical exponent u. As 
is clearly seen from this analysis, the problem belongs to the same universality class of 
random percolation regardless of the value of K considered. 

The same procedure described above can be realized for different values of K as is 
shown in figure 3(6) just for K = 0, 1 and 2.5. Thus, for a given value of K, all the 
curves used in the experiment (for different values of L) collapse into a universal curve 
according to the theoretical prediction. However, R^ is not only a function of 9 and L 
but also of K. As can be seen, the collapsing function is different for each value of K 
considered. This fact determines that the scaling function R^ is not a universal function 
with respect to the variable K (each value of K is represented by using a different type 
of line as indicated) . 

In order to determine the dependence of R^ with K, the main features of the collapsing 
data have to be considered in the range of K between —2 and 3. As can be seen, the curves 
become steeper upon increasing the value of K. In fact, the derivative of the universal 
function R-^ with respect to u behaves as a Gaussian-like function. Thus, we can observe 
that: 

a) the derivatives become more pronounced as K increases. It is possible to establish 
a power law to describe this behavior. Then, 



9^{L) = 9c{oo) + A''L~'/' 



(5) 




(6) 




(7) 
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In a loglog scale the points are very well correlated by a linear function, as in 
equation ([7]), the numerical values of the fitting parameters being very similar for 
the three criteria used here. 

b) the derivatives are narrowed upon increasing K. This behaviour can also be de- 
scribed by a power law according to: 

= CK~^. (8) 

where is the standard deviation of f^^) for each curve. Thus, the standard 
deviation of each derivative versus K when plotted in a loglog scale is very well 
correlated by a linear function (not shown here), with the fitting parameter A = 
4.18 ± 0.02 for A, U and I criteria. 

According to the above equations, a metric factor might be included in the scaling 
function, equation ([S]), in order to collapse all the curves in figure 3(6) onto a single one. 
Following [2B], in figure 3(c) we plot the probability i?^ as a function of the argument 
u' = {9 — 6c) L^^'^K^. As is clearly observed, all the curves collapse onto a single one. It 
is remarkable that more than 2 x 10^ points are included in the collapsing curve. The 
metric factor introduced here, K^, gives an additional indication for the numerical value 
of the exponent A obtained in equation ([H]). 

4 Phase diagram and thermal finite-size scaling 

By using the scheme discussed above, the critical curve, 6c versus K, separating the 
percolating and non-percolating regions, is presented in figure 4. In the studied range, 
three regimes can be distinguished: {€) for K < —1.76, 6c remains constant as K is 
decreased; (ii) from K ^ —1.76 up toK ^ 1.76 {K ^ —1.76 {K ^ 1.76) being the reduced 
critical temperature for the condensation (orderdisorder) phase transition occurring in 
the system), 6c increases linearly with K and {Hi) for K > 1.76, 6c remains constant 
as K is increased. This behaviour can be explained by simple geometrical arguments. 
Namely, lateral attractive interactions favour the nucleation, which in turn increases the 
local connectivity and diminishes the percolating fraction of occupied sites. In contrast, 
repulsive couplings avoid the occupation of nearest neighbour sites, and consequently, 
increase the percolation threshold. In the limit cases, once Kc is reached, the adlayer 
does not vary significantly as \K\ is increased, and 6c reaches its saturation value. Thus, 
6c = 0.518 for K < -1.76 and 6c = 0.662 for K > 1.76. 

As can be seen from the phase diagram, there exists an alternative route for deter- 
mining the critical curve. In fact, the surface coverage can be kept constant while the 
reduced temperature, K, is varied. This procedure does not constitute the standard tech- 
nique for calculating percolation features and it will be revised in detail in what follows. 
It is important to emphasize that the idea of studying phase transitions upon varying a 
controlling parameter and at the same time to keep constant the surface coverage is not 
new. In fact, this strategy has been used in several works related with either kinetic (for 
example, in combination with the well-known ZGB model, etc) [29| [301 13T] or equilibrium 
phase transitions [32j. However, to the best of our knowledge this study has not been 
done for the present model despite its simplicity. 



6 



In figure 5 typical curves of the probabihties Rj^{K) (squares), Rl{K) (circles) and 
Rf^{K) (triangles) are shown for two different values of coverage {9 = 0.5927, full symbols 
and 6 = 0.638, open symbols). For each case (fixed values of 6 and a given criterion), just 
three different lattice sizes are shown in the figure. The curves cross each other in a unique 
point, R^* , which depends on the chosen criterion. These points do not change their 
numerical values regardless of the used coverage. These findings encourage for considering 
a finite-size scaling analysis with the temperature as independent variable by following the 
same rules as in section [3] where the percolation probability was calculated as a function 
of coverage. Thus, for each curve, K^{L) is determined by least mean-square fitting. An 
extrapolation of K^{L) towards the limit L — > cx) by using the scaling hypothesis can be 
argued as 

K^{L) = K,{oo) + C^L-i/^ (9) 

where u is the critical exponent associated with the correlation length and is a 
constant for each criterion used. Equation ([9]) is supported by numerical results as is 
shown in figure 6 for surface coverage 6 = 0.626 as an illustrative example. Thus, all the 
curves (different criteria) are well correlated by a linear function with a quite similar value 
for the ordinate in the limit L ^ oo. It is important to emphasize that K^{L) is almost 
independent of the lattice size. Therefore, the crossing points (figure 5) are located at 
very well-defined values in the A'-axes determining the critical percolation threshold for 
each coverage. 

By using standard finite-size scaling, it is possible to determine the critical exponent 
u which results to be equal to = 4/3 regardless of the value of 9 considered as both (a) 
in the case of random percolation and (b) in section [3l 

As a consequence of the above results, it is possible to collapse all the curves in 
figure 5 onto a single one for each coverage when R^{K) is plotted as a function of 
z = {K — Kc) L^^'^ . As an example, we plot R^ as a function of z for the surface coverage 
9 = 0.638 in figure 7, giving an additional indication for the calculated numerical value 
of the critical exponent u. The same procedure described above was realized for different 
concentrations. This fact demonstrates that the scaling function depends not only on the 
variable K but also on the coverage. 

Thus, for determining the dependence of R^ with K, the main features of the collaps- 
ing data must be considered: the curves become steeper upon increasing the coverage. 
In fact, the derivative of the function R^ with respect to z behaves as a Gaussian-like 
function. Thus, we can observe the following facts: the derivatives (a) become more pro- 
nounced and (b) are narrowed as the coverage increases. The latter is the most important 
finding where one wishes to find the appropriate parametric factors in order to collapse 
all the curves, for different values of 9, into a single universal curve for each criterion. 
Such behaviour can also be described by a power law according to 

Af oc 9-^. (10) 

where is the standard deviation of 

According to the above equations, a metric factor might be included in the scaling 
function in order to collapse all the curves onto a single one for each criterion. In figure 8 
we plot the probability R^{K) as a function of the argument z' = {K — Kc) L^^'^9^. As 
is clearly observed, all the curves collapse onto a single one for each used criterion. This 
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fact allows us to determine the numerical value of the parameter A which results to be 
equal to A = 1.45 ±0.05 regardless of the criterion used. The curves nicely collapse in the 
close vicinity of z' = (close to the critical point) and a tiny deviation is observed as \z'\ 
increases. The scaling analysis given above should be rigorously valid only for sufficiently 
large L and from T in the asymptotic critical regime. However, it can be applied to the 
entire range of L and T if the data fall in the domain of attraction of a simple fixed point 
characterizing only one universality class of critical phenomena. 

It is remarkable that more than 2 x 10^ points are included in the collapsing curve. 
The metric factor introduced here, gives an additional indication for the numerical 
value of the exponent A obtained by using equation ( |T0|) . 



5 Conclusions 

We presented a model to investigate the process of adsorption of interacting monomers 
on a square lattice and studied the percolating properties of the adsorbed phase. By 
using Monte Carlo simulation and finite-size scaling theory, we obtained the percolation 
thresholds for different values of concentration and temperature. From this analysis, 
a critical curve in the 6T space was addressed. The coexistence line, separating the 
percolating and nonpercolating regions, is characterized by three regimes: (1) for K < 
—1.76, 9c remains constant {9c = 0.518) as K is decreased; (2) from K ^ —1.76 up to 
K ^ 1.76, 9c increases almost linearly with K and (3) for K > 1.76, 9c remains constant 
{9c = 0.661) as K is increased. 

Each point in the critical curve was corroborated by following an alternative scheme: 
the surface coverage can be kept constant while the reduced temperature, K, is varied. 
This study does not constitute the standard procedure for calculating percolation prop- 
erties. The results in this paper show that the new technique seems to be a promising 
method for describing the percolation behaviour of adlayer at equilibrium or what we 
called thermal percolation. 

In all considered cases, the finite-size scaling study indicates that the model belongs 
to the universality class of the random percolation model. 
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Figure Captions 



Fig. 1: Fraction of percolating lattices as a function of the surface coverage. Different criteria are 
used for establishing the spanning cluster, namely, R^{9) the probability of finding either 
a rightward or a downward percolating cluster (circles); Rl{p) the probability that we 
find a cluster which percolates both in a rightward and in a downward direction (squares); 
i?f (p) = i \Rf:{p) + R^{p)] =1 IrUp) + Rl{p)] (triangles). Three different values of K 
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are used as it is indicated. Horizontal dashed lines show the i?^* universal points. Vertical 
dashed lines denote the percolation threshold, 9c in the thermodynamic limit L — oo. 
The error bars are smaller than the symbol size. 

Fig. 2: Extrapolation of 9c{k) towards the thermodynamic limit according to the theoretical pre- 
diction given by eq.(l5]). Squares, triangles and circles denote the values of 6c{k) obtained 
by using the criteria /, A and U , respectively. Different values of K are presented as 
indicated. The error bars are smaller than the symbol size. 

Fig. 3: (a) Collapsing plot of the curves for the fraction of percolating samples as a function of u 
For the case K = Each symbol denotes a different value of L (L = 16, 32, 48, 64, 80, 96 
and 112). Each of the solid lines which are simply a guide for the eye represents one of 
the criteria U , A and J, discussed in the text. (6) The same as in (a) for different values 
of K as indicated. The lines are just representative curves for the different criteria, (c) 
The probability i?^ as a function of the argument u' = (9 — 9^) L^^'^K'^, where the metric 
factor i^'^ is included in order to collapse all the curves in Fig. 3a onto a single one for 
each criterion. 

Fig. 4: Phase diagram, 9c vs. K, which shows the curve separating the percolating and not 
percolating regions. Vertical dashed lines at K = —1.76 and K = 1.76 denote the reduced 
critical temperature for the phase transition occurring in the adlayer phase for attractive 
and repulsive interacting particles respectively. Horizontal dashed lines at 9c = 0.518 and 
9c = 0.662 are the critical coverage at saturation regime for K < —1.76 and K > 1.76 
respectively. The error bars are smaller than the symbol size. 

Fig. 5: Typical curves of the probabilities R{^{K) (squares), R^{K) (circles) and R^j^{K) (tri- 
angles) are shown for two different values of coverage. 9 = 0.5927 corresponds to full 
symbols while 9 = 0.638 is denoted by open symbols. For each case (fixed values of 9 and 
a given criterion), just three different lattice sites are shown in the figure (L = 32, 48 and 
64). The universla points R^* are denoted by horizontal dashed lines while the critical 
temperatures in each case are represented by vertical dashed lines. The error bars are 
smaller than the symbol size. 

Fig. 6: Extrapolation of Kq{L) towards the thermodynamic limit according to the theoretical 
prediction given by equation [9] for surface coverage 9 = 0.626. Squares, triangles and 
circles denote the values of Kc{L) obtained by using the criteria /, A and U respectively. 
The numerical value of Kc{oo) is indicated in excellent agreement with the phase diagram 
shown in figure 4. The error bars are smaller than the symbol size. 

Fig. 7: Collapsing plot of the curves for the fraction of percolating samples R^{K) as a function 
of z = {K - Kc)L^''' for the case 9 = 0.638. Each symbol denotes a different value of L 
{L = 16,32,48,64,80,96 and 112). Each solid line, which is simply a guide for the eye, 
represents one of the criteria, U, A and /, discussed in the text. 

Fig. 8: The probability R^ as a funtion of the argument z' = (K — Kc)L^^'^9^, where the metric 
factor 9^ is included in order to collapse all the curves onto a single one for each criterion. 
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